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SUMMARY 


The problem of developing guidance information for changing the orbit of 
a vehicle by using its aerodynamic Lifting capability is considered. Aero- 
dynamic maneuvers reduce propulsive control system requirements to achieve a 
range of orbit transfers of practical interest. As a consequence, rocket 
fuel, and hence, the tremendous cost of transporting additional fuel mass to 
orbit are saved. In order to achieve these savings, guidance laws must be 
developed for the aerodynamic portion of the maneuver sequence. When the 
atmospheric maneuver is accomplished entirely aerodynamically . it is termed an 
aeroglide maneuver. The approach taken for developing the aeroglide guidance 
law is to analytically approximate the solutions to those control problems 
which optimize the final path angles for a given energy loss. A detailed 
analysis of the optimal heading angle problem is provided, but the methods 
used are equally applicable to the coplanar. orbit transfer problem wherein 
the final flight path angle is optimized. 

The optimal heading angle problem for a given energy loss is equivalent 
to the minimum energy loss problem that accomplishes a desired change in the 
orbital inclination of the vehicle. Analytic expressions for the optimal 
controls (bank angle and lift) are developed as solutions to the approximate 
state/Euler system of differential equations for the optimal heading angle 
problem. The optimal control solutions are characterized by three approxima- 
tions, valid in separate regions of the flight, which are derived using asymp- 
totic theory of linear differential equations containing a small parameter. 

The perturbation parameter depends on the scale height of the atmosphere, 
assumed to be exponentially varying with altitude, and the planet's radius. 

The optimal control solutions are a composite of two slowly varying (outer) 
solutions, that are valid in the region of flight near the boundaries, and a 
rapidly varying (inner) solution that is valid in the region where the minimum 
altitude for the flight occurs. Numerical analyses, required to determine the 
matching conditions that continue the two outer solutions through the inner, 
transition zone are not included. Nevertheless, the analytic formulas for the 
optimal controls form the basis for an on-board guidance algorithm and are 
also useful for developing engineering insight into the optimal steering 
policy . 


Aerodynamic heating constraints, neglected in this analysis, need to be 
imposed to obtain useful guidance laws. Additionally, guidance corrections 
for off-nominal atmospheres must be made for on-board use. 
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denotes solution to left of turning point 
denotes solution to right of turning point 
denotes value at max (L/D) 

denotes differentiation with respect to nondimensionai arc length 


INTRODUCTION 

Aerodynamically assisted, orbit transfers are known to be more fuel 
efficient than all propulsive (Hohmann) transfers (ref, 1,), One application 
for aero -assisted maneuvers is a high - Ear th - orbit to low-Earth- orbit transfer 
for returning vehicles, that service satellites, to the Shuttle, Another 
application is for the Mars mission where aero - assisted maneuvers can be used 
for rendezvous operations with other vehicles. NASA soon plans to launch an 
aero-assist flight experiment to make measurements of the environment 
surrounding the vehicle during sustained high altitude, hypersonic flight. 

Recent numerical studies have appeared in the literature to determine 
atmospheric trajectories that will minimize the amount of fuel needed for an 
orbit transfer. In reference 2 optimal atmospheric trajectories are presented 
for a coplanar transfer in which the apogee of the orbit is lowered. In 
reference 3 optimal trajectories are computed for changing the orbital incli- 
nation of a vehicle. Implementation of optimal control policies on-board a 
vehicle, to take advantage of the aero-assist maneuver's fuel efficiency, 
requires either the storage of optimal controls and trajectories for a large 
number of boundary conditions or rapid on-board computations of optimal 
trajectories for specific boundary conditions corresponding to the desired 
mission. Thus, either large storage requirements or fast, on-board computa- 
tional capabilities are necessary to implement optimal maneuvers. An alter- 
native to either of these methods of implementing optimal steering is to 
approximate the optimal control problem so that it may be solved analytically. 
If sufficiently accurate guidance solutions can be found, the demands upon on- 
board computational capability are reduced. 

In references 3, 4, and 5, different approximations have been made to 
determine guidance laws for the orbit plane change problem. In reference 3, 
the term appearing in the differential equation for the flight path angle, 
known as Loh ' s "constant,” is assumed to be constant on an optimal trajectory, 
which makes possible the derivation of an optimal guidance approximation. In 
reference 4, a regular perturbation approach is used to derive optimal guid- 
ance approximations by first neglecting Loh 1 s "constant" and then correcting 
the solution by allowing variation in Loh ' s constant. In reference 3, a 
singular perturbation approximation is derived which simplifies the optimal 
control problem by reducing the order of the state/Euler system, resulting in 
a relatively simple guidance solution. The approach taken herein *is similar 
to those cited above in that the equations of motion -are written in terms of a 
small parameter so that asymptotic theory for differential equations can be 
applied to arrive at a solution to the guidance problem. This research was 
undertaken to improve the accuracy of the approximations cited above. 
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MODELING AND PROBLEM STATEMENT 


The problem of developing a guidance scheme for transferring a vehicle 
from one orbital plane above the Earth to another, using a minimum of fuel, is 
considered. If heating constraints on the vehicle are not imposed, the 
maneuver sequence consists of exoatmospher ic portions, during which the 
vehicle maneuvers using its rocket motor, and an atmospheric portion, during 
which the vehicle maneuvers aerodynamically to change heading by modulating 
angle of attack and bank angle. The imposition of heating constraints on the 
vehicle during atmospheric flight may lead to additional thrusting arcs to 
avoid overheating, reference 6, but such constraints are not considered in the 
present formulation. The maneuver sequence begins with a deorbit impulse to 
cause the vehicle to descend toward the Earth's atmosphere. After entering 
the atmosphere, the vehicle maneuvers to change its heading to the desired 
orbital inclination by modulating angle of attack and bank angle. A boost 
phase follows the aerodynamic maneuver to increase altitude of the vehicle to 
a desired apogee. Finally, a rocket burn is used to circularize the orbit. 
Because the exoatmospheric portions of the flight can be calculated in closed 
form using Kepler's equations and AV impulses to approximate fuel usage, 
interest centers on the atmospheric phase of flight for obtaining analytical 
expressions that can be used for guiding the vehicle. The atmospheric maneu- 
ver problem can be stated in terms of the AV impulses of the minimum fuel 
problem, resulting in a one -dimensional parameterization of the atmospheric 
guidance problem in terms of the deorbit AV impulse and a performance index 
which minimizes the energy lost during the aerodynamic turn, reference 3, The 
vehicle's motion is described in terms of six state variables: specific 

energy, E = V 2 /2 - |i/r; flight path and heading angles, yand down range and 
cross range angles, 9 and <j>; and the radial distance, r of the vehicle from 
the Earth's center. The controls for the vehicle are lift, L; and bank angle, 
jl. The velocity of the vehicle is denoted by V, |i is the universal 
gravitational constant for the Earth, m is the vehicle's mass, and D is the 
drag force on the vehicle, which is assumed to be parabolic with respect to 
lift. Time, t, is replaced by nondimensional arc length, 

s = J (V/r) dt (1) 

as the independent variable. The point-mass equations of motion of the 
vehicle may be written with respect to a vehicle centered reference frame that 
rotates above a spherical Earth, inertial frame as depicted in figure 1. 


dr/ds = rsiny (2) 

d9/ds - cosyc os\|//c os<|> (3) 

d<)>/ds = cosysiny (4) 

dE/ds = -Dr/m (5) 

dy/ds * Lrcos^/(mV 2 ) - [ |i/ ( rV 2 ) - 1 ] cosy (6) 

d\|//ds = Lrsinji/ (mV 2 cosy) - cosycosy tan<[> (7) 


The cross range angle, <]>, is small near a node so that the last term on the 
right-hand side of equation (7) may be neglected if the turn is made in the 
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vicinity of a node. If the boundary conditions on $ and 0 are not prescribed, 
these state variables are ignorable since neither couples into the remaining 
equations. Thus, the order of the system is reduced by two. Since § is 
small, the equation for orbit inclination, 

cosi = cos<(> cosy £ cosy (8) 

and the heading angle approximates the inclination angle near a node. 

r' - rsiny ^ 

E' = - Dr/m ( 10) 


y' = Lrsin(i/(mV 2 c osy) (11) 

- r - , (12) 
y' = Lrcosll/mV 2 - [ |i/(rV 2 ) - 1 ] c osy 

The following nondimensional variables are introduced to determine if any 


further simplifications can be made. Let: 

* 

Z - p r SC L r r e~ h /2m 
h - P(H - H r ) 

where H is the altitude of the vehicle above the Earth's surface, H r is the 
reference altitude for the exponentially varying atmosphere, is the air 
density at H r f S is the vehicle's reference area, and 

r r = R + H r ( i5 ) 

where, R is the radius of the Earth and 

C L = Cl at max (L/D) (16) 

The drag coefficient is assumed to be parabolic 

CD = C Dq + KC 2 l - C Do (l+X 2 ) = C* (1+X 2 )/2E* (17) 

* 

where, E* = (L/D) max and X = Cl/C l . Also let 

v = V -/rr /A 08) 

e = Er r /|i - v 2 / 2 - ( 1 + eh)* 1 09) 

r/r r - (1 + eh) (20) 
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be the nondimensional speed, energy, and radius variables, respectively, where 
“1 ^ 

€ 35 (P^ r ) is a small parameter (of order 10'- 3 for the Earth). Substituting 
equations (13) through (20) into equations (9) through (12), one may rewrite 


the governing equations of motion for the vehicle as 

h = e ’ * ( 1+eh) siny (21) 

e' = -Z r e' h (l+X 2 ) (e + 1 + ehe)/E* (22) 

\\f' * Z r e ( 1+e h) Xsin|i/cosy (23) 

y' - Z r e ‘^Xcos|i ( 1+eh) - { 1 +g [h / (e+1 ) -h] * 1 ) / [2 (e+1 ) ] cosy (24) 

to 0 (g 2 ), where the binomial theorem has been used to expand r/r r . After 
making a final change of variables 

u - In (e+1 ) (25) 

and substituting equation (25) into the system. (21) through (24). one obtains 
to 0 (g) the following approximate model, describing the vehicle's motion: 

u ' = -Z r e- h (l+?i 2 )/E‘ (26) 

Y* “ Z r e ‘^Xsinjl/cosy (27) 

y = Z r e'^X.cosp. - qcosy (28) 

h" = e'^d+eh) siny (29) 


where , 


q = e ~ u /2 - 1. (30) 

The approximate model of the vehicle, (26) through (29). retains the term, q, 
which is a factor in Loh r s "constant" but does not account for variations in 
gravitational acceleration due to altitude. Additionally, the model assumes 
that the effect of altitude changes on the motion are dominated by changes in 
air density rather than changes in potential energy. Thus, potential energy 

changes are neglected on the right-hand sides of equations (26) through (29). 

* 

Finally, variations with Mach number of the aerodynamic coefficients, and 

E* . have been neglected at this order of approximation so that Z r is regarded 
as constant. 


OPTIMAL CONTROL PROBLEM 

If the final path angles are specified, a minimum energy loss trajectory 
to a specified final altitude will result in minimum AV fuel burns to achieve 
a desired final orbit. By restating this problem. in terms of a specified 
energy loss, one may instead, optimize the final path angles. The purpose of 
the restatement is to simplify the final boundary conditions on the optimal 
vertical lift component, which equals zero in this instance. Only the plane 
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change problem is discussed in detail herein. The plane change problem is 
equivalent to extremizing the final heading angle for a specified energy loss, 
so that the payoff is 

J = Vf “ J y' ds. 

The Hamiltonian for this problem is 

H = X 0 y ' + X u u' + Xhh + Xyl . 


(31) 

(32) 


where , 


+ 1 :=> max => max i 

*o-{ <”> 

-1 min \j/ f =>min i 

By choosing X c - + 1 , one seeks the optimal bank angle and lift controls to 
maximize the final heading angle with starting conditions, u e « h^ f Ye * 
final conditions uf , hf , so that the final path angles are optimized. 
Substituting the state equations, (26) through (29), into equation (32), one 
obtains 


H = Z r e "^Xsin|i/c osy - X u Z r e ' ^ ( 1 +X 2 ) /E 
qc osy) 


n- \\ y e ‘ 1 ( 1 + eh) siny + Xy (Z^e “^Xcosji - 

(34) 


Since the final arc length is not specified, H * 0 on an extremal. The optimal 
controls are given by 

3H/3X - -2X U X/E* + sinn + Jlycos|l « 0 (35) 

and 

3H/3|i * cosp. - Xy s in(j. = 0. 
where second-order terms in y have been neglected. 

Solving equations (35) and (36) simultaneously . one obtains 


X = E * V 1 +Xy 2 / ( 2 X u ) with X u > 0 (37) 

cotjil = Xy ( 38 ^ 

The Lagrange multipliers satisfy 

X u — - 3H/ 3u='Xye’ u /2 (39) 

Xy — -3H/dy = - 6 “ 1 Xh(l- f_ eh)cosy- (Xyq + ZXsinpsec 2 y) siny (40) 

X h ' = -3H/3h - H - Xh^ ’ 1 (l+eh)siny+ qX Y cosy - Xh siny 
- -Xh € " ^l+t (h+l)j siny + qXy cosy 
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since H * 0 . Differentiating equation (40). one obtains 

eXy" = -sinyXy - qXy + 0(e) + 0(sin 2 y) (42) 

The solution to equation (42). rewritten as 

e Xy + y Xy + q Xy =0 (43) 

for small y. gives the cotangent of the optimal bank angle control. Equation 
(43) is regarded as a linear equation with coefficients depending on the inde- 
pendent variable, s. The behavior of the solution depends on the coefficient 
of the damping term, which is the path angle, y. According to reference 7. if 
y(s e ) <0. there will be no boundary layer for the solution of Xy at s e . and 
that if y(sf ) > 0. there can be no boundary layer for Xy at Sf. In the 
present problem y(s e ) < 0 for entry into the atmosphere and y(sf) > 0 for 
exit. Thus, no boundary layers are expected in the solution of (43). Since y 
(s G ) < 0 and y (sf) >0. y(st) = 0 for some St» such that s e < St < sf. The 
simplest possible situation is that y vanishes only once on the interval so 
that 


y < 0 for s 0 5 s < st 

y > 0 for st < s < Sf 

y (st ) =0. (44) 

The solution to equation (43) can be sought in two steps with the aid of 
assumptions (44). The first step is to find the solutions for s away from st , 
which break down at s = St- The second step is to find a solution in a neigh- 
borhood of s = St which is asymptotic to the two solutions found in the first 
step. Problems that contain internal transition zones, such as this, are 
known as turning point problems, reference 8. with the turning point, in this 
case, occurring at s = St when y =■ 0, 


ANALYSIS OF THE TURNING POINT PROBLEM 

The optimal bank angle is found from the solution to equation (43) , 
which can be placed in a standard form whose solution is known. For 
convenience later, the independent variable is translated by letting 


x = (s - St) 

so that x(0) * -St and x(st) 351 0. The d 
from Xy to w by letting 

Xy = w exp ( - Jq (y/ 2e ) dt) . 

Substituting equation (46) into equation 
terms, one obtains 


(45) 

pendent variable is transformed 

(46) 

(43) and factoring off exponential 


w" - (26 ) ' 2 Q(x)w = 0 


(47) 
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where . 


Q(x) - y 2 + 2e Y - 4 G q • 


(48) 


Equation (47) is the standard form for analyzing turning point problems with a 

_ 2 

large parameter. (2e) , reference 8. 

Outer Solutions 

The general solution to equation (47) is given by 
w= lQ(x)l [ Aexp (J^(VQ/2€)dx ) + Bexp ( - J* ( Vq/ 2e ) dx ) ] + 0(e) (49) 

where A. B. and c are constants to be determined. The fractional powers of 
Q(x) in (49) are determined using the binomial expansion, assuming Y(x) is 
away from zero, outside the transition zone. From equation (48) 

Vq = y + ey'/y- 2 eq/y (50) 

|Q I 1/4 “ VI Y ! + 0(e ) . (51) 

Substituting equation (46) into equation (49) . together with (50) and (51) , 
one obtains for the zeroth-order approximation to Xy to the left of the 
turning point 


Xy =a 1 E' 1 + a 2 e‘(h-h ) E /y. x <0 

• n 


(52) 


where , 

E = exp (t (q/y) dt ) . 

To the right of the turning point, 

\ = b i F -1 + b 2 e*< h - h mln ) F /Y: x > 0 

where , 

F = exp (if (q/Y)dt) . 


(53) 


(54) 


(55) 


The zeroth-order solutions for Xy given by equations (52) and (54) describe 
the cotangent of the bank angle control outside the transition zone, i.e., for 
lx I > 0. 

The solutions may be regarded as a mixture of terms (aj.E'*, a 2 E/y. b^F’*, 
b 2 F/y) which are functions of the slow scale x and of the term, 
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which is a function of the fast scale, x/e, since h-h m in ~ |(y/e)dx. The 
solutions (52) and (54) break down in the transition zone where y= 0. When x 
is small in the transition zone, y is assumed to have the Taylor series 
expansion 

I x + . . . . = ax + . . . (56) 

1 1 1 x=0 

where, a > 0. Likewise 

y — > a and q ->q(0) h q t (57) 

as x — >0~. Under assumptions (56) and (57), the differential equation for w, 
in (47) , becomes 

W 0 "- (2e)~ 2 [ (ax) 2 + 2e a - 4eq t ]W 0 « 0 (58) 

where, W 0 denotes the zeroth-order inner variable for the solution to equation 
(47) when x is small. Likewise, under the small x assumptions, the outer 
solutions are, from equations (52) and (54): 

X? = a< exp ( * Jx (q r /ax)dx) + a, (ax) * l e - ax? ,2e exp ( J* (q t /ax)dx) (59) 

*o 1 e e 


= bqexp (J* f (q t /ax)dx + b2 (ax) ~ V a:<2/2e exp ( q t /ax) dx ) . (60) 

The inner solution to equation (58) can be expressed in terms of the stretched 
dependent variable, X = Va/e x. This solution leads to an inner expansion 

for denoted by A v . The inner solution. A v , for large X can be matched 

T *o *o 

with the outer solutions, equations (59) and (60), for small x. 


Inner Solution 

Equation (58) for the inner variable W Q can be rewritten in terms of the 
stretched dependent variable X = Va/e x as 

d^W Q /dX^ - (X^/4 + 1/2 - q t /a) W Q ~ 0. (61) 

Equation (61) is the parabolic cylinder equation whose solution is known in 
terms of the parabolic cylinder functions, D-^(X) and D\)(-X), when is not a 
nonnegative integer. The solution to equation (61) is given in reference 9 as 

W 0 (X) = Cl D u (X) + c 2 M-X) (62) 

where. \) = q^/ct-l, so that the solution for Ay 0 is given by 
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Ay o 00 


e 


(63) 


-X 2 / 4 [ Cl D u (X)+ c 2 D v (-X)] . 

The asymptotic properties of the parabolic cylinder functions for large X are 
used to match with the outer solutions for small x. The asymptotic properties 
of the parabolic cylinder functions are given by 

2 

D U (X) - (X) u e' x /4 as X -» + «> (64) 

Dy(-X) - (X)^" 1 e x2/4 'f2n/r{- X) ) as X -> + °° (65) 

where T denotes the gamma function. 


Matching and the Formation of. a Uniform-Expansion f o r Xy (ref. IQ) 

The inner expansion, equation (63). is first expressed in terms of the 
asymptotic limits of the parabolic cylinder functions as X — > ±°°and then 
matched with the outer expansions for small x as x 0-. From equations 
(63) through (65) as X — > + «>, 

A Yo (X) = c^X^t'a-VX 2 / 2 + c 2 (X)- < lt/« VS/r(l-q t /a). (66) 


Since the first term on the right-hand side of (66) is exponentially small 
compared with the second, the first term can be neglected in comparison with 
the second. Next, the inner limit may be expressed in terms of the outer 
dependent variable, x, in preparation for matching, as 


A Yo ~ c 2 C ‘ 


■ q /a 

t 


as X 


(67) 


where , 


C = (VaTe ) q t /a [r(l-q t /a)]/V2^ 


( 68 ) 


Similarly, to the left of the turning point the inner limit is given by 

Ayo~ c i C-l(-*rV«. as X--> - oo. (69) 

r 

Equation (67) is matched with Xy as x -» (r . and equation (69) is matched 

Q r 

with Xy as x — > O', The second term in equation (60) for Xy is, for small 
o 

e , exponentially small compared to the first so that the second term may be 
neglected in comparison with the first. Thus. 

Xy - b 1 e -(q Wnl*/* f l . as x ->0 + (70) 

o 
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Similarly, from equation (59) 


Jf -a, e Aq/a)ln[-x/l-x)] ■ as x ->0\ 
Yo 1 * 


(71) 


A digression to discuss the boundary condition on X is necessary before 

Yo 

r 

proceeding with the matching process. The boundary condition on X^ is zero 

since the final value of y is unspecified. This boundary condition can only 
be satisfied asymptotically for large x (since b^ *0). This is possible 

since q > 0 and y > 0 to the right of the turning point so that the outer 
r 

solution, X , from equation (54) decays from the turning point. Thus. b-. is 
Yo 

r 

simply a scale factor for X y . Thus, the upper limits on the integrals in 

r o 

equation (60) are set equal to + 00 to indicate that the final boundary 

r 

condition on Xy is satisfied asymptotically, i.e., X (x) — ) 0 as x and 

* Yo 

b^ =1. Returning to the matching procedure, one matches equation (54) with 
equation (67), using (70), by equating the two .limits as X — > + and x — » 0 " h : 


C 2 -> C exp (q/y~ q t /otx)dt. 


(72) 


so that 


c 2 = C exp 


ft - 00 


(q/y- q t /at)dx. 


(73) 


Similarly, to the left of the turning point, equation (52) is matched with 
equation (69), using (71), by equating the two limits as X and x -> 0 : 


c 1 


ai ( ~x e ) Cexp 



q /ax) dT . 


(74) 


so that 


c x = a x (-x e ) qt/ “cexp £ e ( q /y - q t /or) dt . (75) 

Equations (73) and (75) provide the connections between the coefficients, ci 
and C2. of the inner solution and the coefficients of the outer solutions, ai 
and bi, that are determined by the boundary conditions. Thus, the two outer 
solutions have been continued through the transition zone. The matching 
procedure showed that exponentially small terms were neglected in forming a 
leading order uniform approximation. Thus, a 2 = b 2 ^ 0 for the leading order 
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approximation. These terms enter into the solution at the next higher approx- 
imation (ref. 10). The leading order uniform approximation for is con- 
structed by combining equation (66) with the connection formulas. (72) and 
(74): 


Uniform Approximation for 

2 

Xy u = C e~ ax ( A(x)Dx>(Va/ex) + B (x) D u ( Va/e x) ) 

where , 


A (x) * a jL (-x e )\ /a exp J* e (q/y- q t /«T)dl 

(77) 

f+ oo 

B (x) * exp J x (q/y- q /at) dr 

(78) 

C =(Va7€) q t /a [r(l-q t /a)] /V2jT 

(79) 

u = q t /a - 1 = - (n t - 2 ) / (nt - 1) 

(80) 

nt “ (L cos|i/(W-mV^/r r ) x = o- 

(81) 


In equation (81), n t is the "reduced" vertical load factor, evaluated at the 
turning point! that is, n^_ is the ratio of the vertical lift force at the 
turning point to the weight of the vehicle, reduced by the centrifugal force 
at the turning point. 

The results of this section are summarized by equations (76) through 
(81) and the following equations for the outer and inner solutions: 


Outer Solutions: 

r r+oo 

Xy = exp J (q/y)dl, x >> 0 (82) 

o ^ 


e 

Xy =Xy(x G ) exp (- 


J ( q/y) dx) . x << 0 (83) 


Inner Solution : 


'Yo 


= e ax /4e j c (Va/ex) + c 2 D\>( Ww/e x) 


(84) 
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C L = C A (0) 


(85) 


c 2 = C B (0) (86) 

Equations (76) through (81) give the leading order uniform approximation for 
the cotangent of the optimal bank angle control over the entire interval, 
x 0 <x £ +oo. In equation (79), C is a constant for a fixed atmosphere, 

"reduced" load factor at the turning point and desired minimum altitude for 
the maneuver. The coefficients A(x) and B(x) are functions of the slow depen- 
dent variable, x. while the parabolic cylinder functions and the exponential 

factor are functions of the fast dependent variable, X = Va/ e x , in the 
transition zone. The uniform expansion, for large |X| , is asymptotic to the 
outer solutions, which are functions of the slow variable, x, and satisfy the 
boundary conditions on Xy at the two end points. The inner solution, equation 
(84) , depends on the minimum altitude for the maneuver and the "reduced" load 
factor at the turning point. With the uniform approximation for the bank 
angle control in hand, the reader’s attention is focused on determining the 
optimal, normalized lift, X. 


Determination of Normalized Lift 

The approximation for the multiplier, Xy , has been derived in terms of 

/* 

the state and dependent variables, so that Xy can be determined, also, in 
terms of the state and dependent variables. The Hamiltonian can thus be 
expressed in terms of the state and dependent variables and the unknown 
multiplier, X u , needed to determine the optimal normalized lift, X. 

The Hamiltonian can be rewritten in terms of Xy and (2X U /E*) as 



using equations (34), (37), (38), and. (40). In equation (87) 

b = Z~ 1 (y\y u + X Yu q) (88) 

Z - p r SC L r r e' h /2m (eq. 13) (89) 

^ = Ce' ax2/4e {A(x)D u (Va/e x) + B(x)D u ( Wa/e x) ) (eq.80) (90) 

Thus, (2X U /E*) is determined as the positive root of equation (87), 

(2X u /E*) = -b + yfb 2 + (1 + Xy u 2 ) (91) 
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Thus, the approximation for the optimal normalized lift is. from equation 
(37). 


X = a + "V a 2 + 1 (92) 

where , 

a = b/yj 1 + Xy u 2 (93) 

The variable, b, in equation (93) depends on from equation (88). 

Differentiating equation (90) , one obtains 

X'y^ C e { - ( q/y-q^/ax) A (x) D^(X) - ( q /y- q^_/ax) B(x) D^( - X) 

+ A (x) d [e 'ttx 2 /4e D y(X)]/dx + B(x) d[ e -°« 2/4e D^(-X)]/dx (94) 

The derivatives in (94) can be found from the recurrence relation (refs. 11 
and 12) 


d[e'^ 2/4 D^UUl/dS =-e‘^ 2/4 Dx)(U. (95) 

so that equation (94) may be simplified as 

Xy u = - (q/y-q t /ax) Xy u + Ce‘ ax2/4e Va/e[-A(x) Dq t / a (X) + B(x) Dq t / a (-X)] 

(96) 


From equation (88) 

b - Z' l y{ (q t /ax)Xy u + C e' ax2/4e Va 71 f-A(x)Dq t / a (X) + B (x) D /a ( -X) ] ) 

(97) 

so that the coefficient of (2^ U /E*) in equation (87) can be determined, and 
hence, the value of X can be determined from equations (92) and (93). 

The value of X for large |X| can be determined from the asymptotic 
properties of Dq^_/(x(X) : 

D q t / a ( X) -(X)-9 t / «- 1 e x2 /4 /T( - q t /a) 

- Hi -q t /a)Dq t / a - 1 (X) /T( -q t /a) 

“ ( -q t /a)Dq t / a - l (X) , as X -> + » (98) 

Thus, for large positive X, equation (97) has the asymptotic property, 

b -Z' 1 y{(q t /ax)A.y -(q t /ax)A.y u ) =0. (99) 

where exponentially small terms have been neglected. Thus, outside the 
transition zone, b -0, so that 
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X~ 1 


as X —> + «> 


( 100 ) 


which implies 

C L ~ C L* as x + (101) 

Similarly , 

Cl ~C L * as X ->-00 (102) 


also . 


Formulas have been derived for an approximation to the optimal 
normalized lift coefficient. X, by solving a quadratic algebraic equation for 
the unknown multiplier. 2^ U /E* , from which X is obtained from the extremal 
control formula in terms of 2^ U /E* . The value of X was shown to be asymptotic 
to one. outside the transition zone. This implies that the lift coefficient 
is asymptotic to its value at max (L/D) , outside the transition zone. 

The next section contains a summary of the formulas needed to synthesize 
the optimal guidance commands. 


Summ a ry of U plimal Guidance Commands 

Formulas have been derived in the previous two sections for the optimal 
bank angle and lift controls for the problem of maximizing heading angle for a 
given energy loss. The formulas assume that full, reduced - order state and 
dependant variable estimates are available from the vehicle* s navigation 
computer. For preliminary simulation studies one may assume that the reduced- 
order model, derived herein, supplies these estimates perfectly (i.e., 
uncorrupted by noise and other uncertainties) . 

The following is a summary of the optimal control formulas: 


|t=arccot (Xy^) 


where . 


Xy u = Ce' ax2/4e [A(x) Du(Va/e x) + B(x) D\)( Wa/e x )] 
A(x) = Xy(x e ) ( - x e ) < l t / a exp J x e ( q /y- q t /at) dT 


r 


B (x ) « exp J (q/y- q t /ax) dt 


C = (VoTe) V a tr(l-q t /a)l t^2n 


(103) 

(104) 

(105) 

(106) 
(107) 
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( 108 ) 


\)= q t /a- 1 = - (nt - 2)/(nt - 1) 

n t = [L cosM^CW - mV2/r r )] X=Q (109) 

e = l/(Pr r ) (110) 

Dy(X) . D\) ( - X) are given in tabular form (ref. 13) 

x = st + J 0 (V/r)dt : x e = -St (111) 

q = e " u /2 -1 = ( M>V 2 r r - 1) : q t = q(0) (112) 

Y= JVdx (113) 

y « Z Xcosji.~q : ft - Y I x= q (114) 

Z = Z r e' h ; h = p(H-H r ) (115) 

Z r = p r SCL*r r /2m (116) 

Lift Control 

X= a + Va 2 + 1 (117) 

where* 

a = b /J 1 + X Y 2 (118) 

b=Z' 1 Y((q t /ax)Xy u +Ce ax2/4e Va/e [-A(x) D qt / a (Va/e x) + B(x)D q ^ /a ( -Va/e x)] 

(119) 

Values of the State and EfilatfiiLJTaxlahl&a 

h = Jli' dx . u - Ju'd x, \j/=Jy'dx (120) 

h=e'i(l+eh)Y (121) 

u' = - Z (l+X 2 ) /E* (122) 

y'=ZX.sin(l (123) 

E = p(e u -l)/r r (124) 

V = V 2 [E + p7(r r + H)] (125) 
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In order to implement the feedback formulas for the optimal guidance 
commands, one must write a numerical simulation to determine the integrals for 
the state equations as well as the Integrals needed for matching. While this 
has not been completed, herein, an analytical example of the procedure for 
conducting numerical studies is provided in the next section. 


ANALYTICAL EXAMPLE AND DISCUSSION 

Although numerical simulations are necessary to complete the matching 
procedure and to determine the accuracy of the optimal guidance approximation, 
the behavior of the bank angle control may be examined in the transition zone 
where the guidance law reduces to an analytical formula. Data from reference 
3, listed in Table 1. are used for the vehicle model in this example and tabu- 
lated values of the parabolic cylinder and gamma functions, given in reference 
13, are used to calculate the guidance commands. 

Table 2 contains estimates of the guidance parameters as a function of 
minimum altitude for the maneuver. The constants, ci = -.0542 and C2 ** ,2748, 
were selected to form the inner approximation, Ay^ , as a linear combination 

of parabolic cylinder functions multiplied by The index, X) =* -.5, for 

the parabolic cylinder functions was selected for this example, which corre- 
sponds approximately to a 150,000 foot altitude trajectory with a bank angle 
at the turning point of 75° and a normalized lift coefficient, 1. 

The index for the parabolic cylinder functions, depends on the 
vertical, "reduced" load factor, nt* at the turning point. n* is a function 
of H m ^ n , V, and Xcosp. at the turning point, so that care must be exercised in 
estimating \). The values for \) given in Table 2 lie in the range, -1< l)<0, 
which corresponds to a range of reduced load factors, n t > 2. For 1 < n t ^2, 
u is non - negative , and when x> is a non-negative integer, the guidance formulas 
are not valid sine e D-q (X)and D^(~X) are linearly dependent in this case. 
Although the guidance solutions can be formulated more generally in terms of 
other linearly independent solutions of the parabolic cylinder equation, one 
would expect, for most applications, that the reduced load factor exceeds a 
value of two. In this event U is a negative number greater than -1. 

Data for Ay^ and the corresponding bank angle. M 0 , are tabulated in 

0 

Table 3 and plotted in figure 2. Values for the two outer solutions, X y and 

r 0 

r 

X , are tabulated in table 4 as functions of the inner variable, X. The 
*o 

outer solution approximations for small x: 
r 

n 0 - cot' 1 [C2 (X/2) ■ 1/2 ] ; X > 0 (126) 


c 


cot - [ci C - X/2) *1/2] 


X < 0 


(127) 


are also plotted in figure 2 so that the behavior of the inner and outer solu- 
tions can be observed in the overlapping regions, where these solutions match. 
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For purposes of comparison the optimal hank angle control, given in 
reference 3, is superimposed on figure 2. The agreement between the inner 
solution guidance approximation. M 0 . and the optimal bank angle control. |l 0 pt* 
is excellent in the neighborhood of the turning point. X = 0 . The agreement 
between the optimal bank angle, |iopt. and its approximation, ax. near the 
turning point is also good. Away from the turning point the outer solution 

r 2 

approximations for small |x| , [i Q and \i Q , are asymptotic the inner solution M 0 . 

for large |X|. The outer solutions for large X are not in numerical agreement 
with the optimal solution for the bank angle, indicating the need to adjust 

2 

the initial condition on |i 0 by adjusting Xy (x e ) . Matching the inner and 

f o 

outer solutions for specified initial conditions requires a numerical 
simulation, which is not included in this report. 

For numerical studies one could begin by integrating trajectories from 
the turning point to the boundaries using the estimates for the state given in 
Table 2. The integration to the right of the turning point could start at 
x » 0 + . by extrapolating the state from x - 0 to x =» 0 + using a Taylor expan- 
sion, Using the guidance formulas as feedback controls, the reduced state 
equations would be Integrated to determine the change in the state from the 
turning point. The flight path angle must remain positive to the right in 

r 

order to satisfy the boundary condition on Xy At the terminal point the 

*o* 

energy lost and heading angle change from the turning point would be 
available. These values would be added to the values of "delta” heading angle 
and energy loss for a similar state - equation integration to the left of the 
turning point. The integration to the left, however, will undoubtedly not 
match the initial conditions on the state. In this case the initial condition 

2 

on can be adjusted to change the outer solution to the left of the turning 
*o 

2 

point. An interative process ensues because X (x e ) is proportional to c of 

*o 1 

the inner expansion, Ay 0< The numerical analysis would proceed by adjusting 
not only ci but also some of the parameters in Table 2. New values for the 
index of the parabolic cylinder functions and a new estimate for Ay^ would be 

calculated. This procedure would be repeated until the initial conditions are 

r 

satisfied sufficiently accurately and the final condition, X ~ 0 as x — > + 00 

i 0 

with y>0 , holds. A numerical procedure such as that outlined above should 
produce bank angle commands that are qualitatively similar to those computed 
with a numerical optimization computer code. 

To complete the numerical analysis, one would use the uniform guidance 
law, obtained by matching with the reduced model, as a feedback law for the 
full system, equations (2) through (7), to compare optimally guided trajec- 
tories with trajectories computed using a numerical optimization procedure. 

This analytical example has shown that optimal trajectories for the 
minimum fuel/plane change problem can be parameterized by the minimum altitude 
for the maneuver in the following sense. Each minimum altitude corresponds to 
a trajectory that satisfies initial conditions and results in a loss of energy 
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and change in heading angle to the left of the turning point. The energy loss 
and "reduced" vertical load factor at the turning point determine the index of 
the parabolic cylinder functions. The loss of energy and change in heading 
angle to the right of the turning point, when added to the corresponding 
values to the left of the turning point, result in an energy loss and maximum 
(minimum) heading angle change for the entire maneuver, parameterized by the 
minimum altitude of the maneuver. 


CONCLUDING REMARKS 

Formulas have been derived for computing approximations to extremal 
controls, normalized lift and bank angle, for the problem of minimizing energy 
loss to change the orbital inclination of a vehicle. Asymptotic analysis was 
used to derive formulas for the bank angle control that were valid in three 
different regions of the trajectory. Two outer solutions were found to 
approximate the bank angle control near the upper limits of the Earth's atmo- 
sphere. These solutions break down at the minimum altitude of the maneuver, 
when the flight path angle is zero. In the transition zone, near Y“ 0, the 
optimal bank angle control was described in terms of parabolic cylinder 
functions. The two outer solutions were connected through the transition zone 
by forming a uniform solution, in terms of the parabolic cylinder functions, 
to approximate the optimal bank angle control over the entire trajectory. The 
approximate, optimal lift coefficient was subsequently determined by solving a 
quadratic algebraic equation. The lift coefficient was determined to be 
asymptotic to its value at max (L/D) when the trajectory is far from the 
turning point. The guidance formulas are in terms of state and dependent 
variables and parameters that define the index and argument of the parabolic 
cylinder function. A numerical example was given to illustrate the utility of 
the method used to construct the guidance formulas. Construction of the inner 
solution showed that optimal trajectories will be parameterized by the minimum 
altitude selected for the flight and the "reduced” vertical load factor at the 
turning point, which determines the index of the parabolic cylinder function. 
Numerical studies, however, must be conducted to obtain matching conditions 
and to verify the accuracy of the guidance law. The guidance formulas require 
feedback of energy, altitude, and flight path angle from the vehicle's 
navigation computer and are suitable for forming the basis for an on-board 
optimal guidance system. 
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TAB LE JL 




VEHICLE/MODEL 

PARAMETERS 




s - 

126 ft 2 


H=1 

.409 x 10 16 

f -sec " ^ 

m - 

332 slugs 


6 = 1. 

29 10' 3 


V 

^ .151 


Hr = 

180 » 000 ft 


p = 

(27000 ft) ' l 


Pr = 

.356 x 10' 4 

lbm/ f t ^ 

r r 

= 20.900,000 ft 


8c = 

32.2 Ibra/sl 

ug 

Zr 

= C L ‘p r Sr r /2mg c - 

.662 (eq. 116) 

E* = 

2.36 



TABLE 2 


GUIDANCE SYSTEM PARAMETER ESTIMATES 
(Ref. 3) 


^min (Est) 

Z (Hmin) 

V t (Est) 

^t 

(ftxlOOO) 

(eq 115) 

(ftxlOOO) 

(eq 112) 

175 

. 797 

25 

.0787 

150 

2.01 

24 

. 170 

125 

5 . 08 

23 

.274 


n t (4 t “ 

75°) V 

a 

Va/€ 

(eq 109) 

(eq 108) 

(eq 114) 


2.62 

- .383 

.128 

10.0 

3.06 

- .515 

. 350 

16.5 

4.80 

- .737 

1.04 

28.5 
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TABLE 3 


Inner Solution 


Ay Q (X) - [-.0542 D- 1/2 (X) + 

.2748 D- 1 / 2 (-X)]e' x2/4 


M_ (X) = cot ' U Ay ) 
0 Jq 



Multipliers 

Bank Angles 

(deg) 

X Ay ( X ) Ay(-X) 

M(X) 

M( -X) 


o 

o 

.2688 

. 2688 

75.0 

75.0 

0.5 

.3420 

.1620 

71 . 1 

80.8 

1.0 

. 3640 

.0621 

70.0 

86.4 

1.5 

.3433 

- .0050 

71 , 1 

90.3 

2.0 

. 3045 

- .0368 

73.1 

92.1 

2.5 

. 2668 

- .0459 

75 . 1 

92.6 

3.0 

.2375 

- .0455 

76.6 

92.6 

3.5 

.2159 

- .0425 

77.8 

92.4 

4.0 

.1998 

- .0396 

78.7 

92.3 

4.5 

. 1871 

- .0371 

79.4 

92.1 

5.0 

, 1767 

- .0350 

80.0 

92.0 
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XAKUL.A 


Outer Solutions 



Xy = 
To 

Xy = 
*0 

.2748 (X/2)' 1/2 ; 
-.0542 (-X/2)‘ 1/2 ; 

mJ 

^0 

= cot “ * (X^) 
“ cot " * (A^) 



Multipliers 



Bank Angles 

(deg) 

X 

Xy 
1 0 

4 

'O 


u r 

4 

. 25 

.7773 

- . 1533 


52.1 

98.7 

.50 

. 5496 

- . 1084 


61 . 2 

96.2 

1.0 

.3886 

- .0767 


68.8 

94.4 

2.0 

.2748 

- . 0542 


74.6 

93.1 

3.0 

.2248 

-.0443 


77.4 

92.5 

4.0 

. 1943 

- .0383 


79.0 

92 . 2 

5.0 

.1738 

- . 0343 


80.1 

92.0 
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BANK ANGLE (DEG.) 



C c )^n') 1 
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